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In this series of eight papers we present the appHcations 
of methods from wavelet analysis to polynomial approxi- 
mations for a number of accelerator physics problems. In 
this part, assuming a sinusoidal field variation, we consider 
the analytical treatment of the effects of insertion devices 
on beam dynamics. We investigate via wavelet approach 
a dynamical model which has polynomial nonlinearities 
and variable coefficients. We construct the corresponding 
wavelet representation. As examples we consider wigglers 
and undulator magnets. We consider the further modifica- 
tion of our variational approach which may be applied in 
each scale. 

1 INTRODUCTION 

This is the third part of our eight presentations in which 
we consider applications of methods from wavelet anal- 
ysis to nonlinear accelerator physics problems. This is a 
continuation of our results from [l]-[8], which is based on 
our approach to investigation of nonlinear problems - gen- 
eral, with additional structures (Hamiltonian, symplectic or 
quasicomplex), chaotic, quasiclassical, quantum, which are 
considered in the framework of local (nonlinear) Fourier 
analysis, or wavelet analysis. Wavelet analysis is a rela- 
tively novel set of mathematical methods, which gives us a 
possibility to work with well-localized bases in functional 
spaces and with the general type of operators (differential, 
integral, pseudodifferential) in such bases. In this part we 
consider effects of insertion devices (section 2) on beam 
dynamics. In section 3 we consider generalization of our 
variational approach for the case of variable coefficients. In 
section 4 we consider more powerful variational approach 
which is based on ideas of para-products and approxima- 
tion for multiresolution approach, which gives us possibil- 
ity for computations in each scale separately. 

2 EFFECTS OF INSERTION DEVICES 
ON BEAM DYNAMICS 

Assuming a sinusoidal field variation, we may consider ac- 
cording to [9] the analytical treatment of the effects of in- 
sertion devices on beam dynamics. One of the major detri- 
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mental aspects of the installation of insertion devices is 
the resulting reduction of dynamic aperture. Introduction 
of non-linearities leads to enhancement of the amplitude- 
dependent tune shifts and distortion of phase space. The 
nonlinear fields will produce significant effects at large be- 
tatron amplitudes. The components of the insertion device 
magnetic field used for the derivation of equations of mo- 
tion are as follows: 
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with k'^ + ky ^ k^ = (27r/A)-^, where A is the period 
length of the insertion device, Bq is its magnetic field, p is 
the radius of the curvature in the field Bq- After a canonical 
transformation to change to betatron variables, the Hamil- 
tonian is averaged over the period of the insertion device 
and hyperbolic functions are expanded to the fourth order 
in X and y (or arbitrary order). Then we have the following 
Hamiltonian: 
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We have in this case also nonlinear (polynomial with de- 
gree 3) dynamical system with variable (periodic) coeffi- 
cients. As related cases we may consider wiggler and un- 
dulator magnets. We have in horizontal x — s plane the 
following equations 

X = -s—Bzis) (3) 
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where magnetic field has periodic dependence on s and hy- 
perbolic on z. 

3 VARIABLE COEFFICIENTS 

In the case when we have situation when our problem is de- 
scribed by a system of nonlinear (polynomial)differential 



equations, we need to consider extension of our previous 
approach which can take into account any type of variable 
coefficients (periodic, regular or singular). We can pro- 
duce such approach if we add in our construction additional 
refinement equation, which encoded all information about 
variable coefficients [10]. According to our variational ap- 
proach we need to compute integrals of the form 

/ b,,{t){^,f'{2"H-h){^2fH2'"t-k2)dx, (4) 
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where now bij {t) are arbitrary functions of time, where trial 
functions tpi, if2 satisfy a refinement equations: 

= X! "■ik^^i'^t ^ ^) (5) 
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If we consider all computations in the class of compactly 
supported wavelets then only a finite number of coefficients 
do not vanish. To approximate the non-constant coeffi- 
cients, we need choose a different refinable function tfj, 
along with some local approximation scheme 

{Btf){x) Yl FiADvsi'ih - k), (6) 

where Fg ^ are suitable functionals supported in a small 
neighborhood of and then replace 6^ in (Q) by 

B£bij{t). In particular case one can take a characteristic 
function and can thus approximate non-smooth coefficients 
locally. To guarantee sufficient accuracy of the resulting 
approximation to (Q) it is important to have the flexibility 
of choosing (^3 different from (pi, tp2- In the case when D 
is some domain, we can write 

b^3{t) b= E b,jit)xD{2h-k), (7) 

0<k<2' 

where xd is characteristic function of D. So, if we take 
^4 = Xd, which is again a refinable function, then the 
problem of computation of (^ is reduced to the problem of 
calculation of integral 

H{kuk2,k3,ki) = Hik) = 
f ip4{2H - k,)ip3{2h - k2) X (8) 

ipf^(2''t-k3)ipi^{2H-k4)dx 

The key point is that these integrals also satisfy some sort 
of refinement equation: 

2-\^\H{k)^J2^^k-eH{i), Ai = di+d2. (9) 
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This equation can be interpreted as the problem of com- 
puting an eigenvector. Thus, we reduced the problem of 
extension of our method to the case of variable coefficients 
to the same standard algebraical problem as in the preced- 
ing sections. 



So, the general scheme is the same one and we have only 
one more additional linear algebraic problem by which we 
in the same way can parameterize the solutions of corre- 
sponding problem. As example we demonstrate on Fig. 1 a 
simple model of (local) intersection and the corresponding 
multiresolution representation (MRA). 




Figure 1 : Simple insertion. 




Figure 2: MRA representations. 

4 EVALUATION OF NONLINEARITIES 
SCALE BY SCALE 

We consider scheme of modification of our variational ap- 
proach in the case when we consider different scales sep- 
arately. For this reason we need to compute errors of ap- 
proximations. The main problems come of course from 
nonlinear terms. We follow the approach from [11]. 

Let Pj be projection operators on the subspaces Vj,j S 

Z: 

Pj : L^(R) ^ Vj (10) 

(P,/)(x) = E < /, V^j,k > fj.kix) 

k 

and Qj are projection operators on the subspaces Wf 

Q,=P,-,-P, (11) 

So, for u G i^(R) we have uj — Pju and uj £ Vj, 
where {Vj},j G Z is a multiresolution analysis of L^(R). 



It is obviously that we can represent Uq in the following 
form: 

n n 

ul = 2Y,{P,u){Q,u) + Y.{Q,u){Q,u) + ul (12) 

In this formula there is no interaction between different 
scales. We may consider each term of ( p^ as a bilinear 
mappings: 

M'vw ■■ X ^ L\K) = V,®,,>,Wr (13) 



Mi^w ■■ X Wj ^ l2(R) = Wj. (14) 

For numerical purposes we need formula ( p^ with a finite 
number of scales, but when we consider limits j cx) we 
have 

u2 ^ ^{2PjU + Qju){Qju), (15) 



which is para-product of Bony, Coifman and Meyer. 

Now we need to expand (|l2|) into the wavelet bases. To 
expand each term in ( p^ into wavelet basis, we need to 
consider the integrals of the products of the basis functions, 
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where j' > j and 

iplix) = 2~^/^xJj{2~^x ~ k) 



are the basis functions. If we consider compactly supported 
wavelets then 

Ml^'^^(k,k'J) = Q for |fc-fc'|>fco, (18) 

where feg depends on the overlap of the supports of the basis 
functions and 
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Let us define jo as the distance between scales such that for 
a given e all the coefficients in ( p^ with labels r — j — j', 
r > jo have absolute values less than e. For the purposes 
of computing with accuracy e we replace the mappings in 
(0), &>hy 

M^w ■■ yj X -> V, ®j<r<jo Wr (20) 



Since 
and 



Mi^w : X W, ^ V, ®j<r<oo (21) 
yj ®3<3'<3o = V,,-x (22) 
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we may consider bilinear mappings (^0|), (|T]) on V^q-i x 
Vja-i - For the evaluation of (pO|), ( ^l| ) as mappings V^q _ i x 



V^g-i ^ K?o-i we need significantly fewer coefficients 
than for mappings (|2C|), (plj). It is enough to consider only 
coefficients 

f-OO 
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(24) 

where (p(x) is scale function. Also we have 

M(fc, /c', i) = 2-^I^Mo{k ~£,k' - (), (25) 

where 

Mo{p,q) = j ip{x-p)ip{x - q)ip{x)dx (26) 

Now as in section (3) we may derive and solve a system of 
linear equations to find Mo {p, q) and obtain explicit repre- 
sentation for solution. 
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